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We use a power grid model with M generators and N consumption units to optimize the grid and its con- 
trol. Each consumer demand is drawn from a predefined finite-size-support distribution, thus simulating the 
instantaneous load fluctuations. Each generator has a maximum power capability. A generator is not over- 
loaded if the sum of the loads of consumers connected to a generator does not exceed its maximum production. 
In the standard grid each consumer is connected only to its designated generator, while we consider a more 
general organization of the grid allowing each consumer to select one generator depending on the load from 
a pre-defined consumer-dependent and sufficiently small set of generators which can all serve the load. The 
model grid is interconnected in a graph with loops, drawn from an ensemble of random bipartite graphs, while 
each allowed configuration of loaded links represent a set of graph covering trees. Losses, the reactive char- 
acter of the grid and the transmission-level connections between generators (and many other details relevant to 
realistic power grid) are ignored in this proof-of-principles study. We focus on the asymptotic limit, N — > °° 
and N/M — > D = 0(1) > 1, and we show that the interconnects allow significant expansion of the parameter 
domains for which the probability of a generator overload is asymptotically zero. Our construction explores the 
formal relation between the problem of grid optimization and the modern theory of sparse graphical models. We 
also design heuristic algorithms that achieve the asymptotically optimal selection of loaded links. We conclude 
discussing the ability of this approach to include other effects, such as a more realistic modeling of the power 
grid and related optimization and control algorithms. 

I. INTRODUCTION 

The existing power grid is complex and far from being optimized. The anticipated installation of small-scale distributed 
generators and storage devices, as well as the addition of many ancillary backup lines and control devices at both transmission 
and distribution levels, imply that the main future challenges will require intelligent planning, optimization and control of this 
ever growing grid. This optimized and efficient grid of the future, that incorporates new hardware and concepts such as renewable 
generation and distributed storage, has been labeled "smart grid" (DEI- Accounting for many important details of the power 
distribution and transmission is not feasible without much simplification. Simplified models must identify the significant effects, 
extracting and analyzing and later probing each of them separately and in combinations. (See, e.g. related review articles 
01113.) 

We adopt this "discovery though simplification" approach and focus in this paper on improving the functioning and control of 
the electric grid on the power distribution level. Specifically, our prime focus is on preventing overloads of the power generation 
units caused by fluctuations in the demand by efficiently utilizing ancillary lines. This approach is justified in the context of a 
city-scale transmission system, in which many significant loads and generators are in geometric proximity of each other and the 
cost of building new ancillary lines is not prohibitive. We assume that shedding excess load is not an option and we consider 
the possibility of redistributing the load via a system of interconnects over a larger grid than what is used under normal (no 
overload) conditions. Fundamentally, we ask the following question: can an intelligent arrangement of 0(M) ancillary lines 
among a system of M generators, each serving D customers, possibly reduce the probability of a generator failure from a finite 
number to zero in the limit of M — > °° and D= 0(1)7 In fact, an assumption that a drastic improvement is possible stems from 
basic information-theoretic intuition: any finite error probability can be reduced to zero via properly introduced redundancy (6|. 
We show that adding interconnects to the power system is related to adding redundancy to the information system, in the sense 
that achieving the asymptotic overload-free distribution is indeed a possibility for an idealized grid. 

The model of power line considered in the paper mimics how the ancillary power lines are operated on the distribution level 
power grid. A city-scale power grid has an intervene loopy structure, however these loops are typically used for ancillary (backup 
and/or maintenance) purposes and operators aim at avoiding running current over loops. We adopt this strategy and assume that 
the system contain (or may be built with) many switches and that it is operated in the manner that for each given configuration 
of the pre-installed on-off switches the currents flow over trees, i.e. subgraphs of the full loopy distribution graph without loops. 
We also assume that the lines are sufficiently short (kilometers, not hundreds or thousand of kilometers), and thus thermal losses 
are not important and reactive parts of line impedances can be ignored. These assumptions correspond to the so-called DC 
approximation [7|. Combination of the loop-free structure of power flows (for any given configuration of switches) with DC 
approximation allows to model electricity delivery as an abstract commodity flow JS). For completeness of our description, let us 
also note that some other effects of realistic power flows are ignored in our first publication on this emerging subject. Thus, we 
do not consider back flows (consumers turned into distributed generators) and for that matter we do not discuss at all the entire 
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scale up (transmission level) structure of the grid, inhomogeneities and spatio-temporal correlations in loads and generation, 
and all economy, pricing and regulation effects [9|. In essence, our main task here consists in establishing the existence of 
fundamental limits, and bounds on the idealized asymptotic failure-free regime, and followed by developing efficient and simple 
algorithm controlling switching in the grid. 

To solve our model we use the cavity method and its computational realization via population dynamics, introduced in the 
statistical physics of disordered systems ifTOlfTTIl and recently adapted to the analysis of Shannon (phase) transitions in constraint 
satisfaction |12| and error-correction |[T3l [T4l . This method explores the famous fact that the Bethe-Peierls |[T5l [T^l or Belief 
Propagation ifTTl [T8l (BP) scheme exactly solves probabilistic models on graphs without loops (a tree). The method allows to 
evaluate the ensemble average over configurations of allowed discrete switchings on the grid balancing load with generation. 
We study how the number of allowed solutions scale with the size of the grid, and identify a transition (in the space of model 
parameters) from a satisfiable (SAT) domain, in which the number of solutions is exponentially large (in the system size) and 
where BP algorithm finds a valid solution easily, to an unsatisfiable (UNSAT) domain in which the number of valid solutions is 
with high probability zero (thus load shedding would be required). 

Assuming that the two ways communications between consumers and generators exists we design a stochastic local search 
algorithm, coined WalkGrid, that is able to find optimal configuration of switches almost anywhere in the SAT phase. We also 
utilize aforementioned BP analysis and develop a BP-based message passing scheme for efficient search of SAT configuration 
of switches. Generally, WalkGrid outperforms the BP-based algorithm almost anywhere in the SAT phase. Note however that 
one useful feature of BP (not readily available in WalkGrid) consists in its ability to count number of available SAT solutions, 
and thus have a direct algorithmic test of the distance to failure, i.e. distance to the SAT-UNSAT transition. 

The material in the manuscript is organized as follows. We introduce our model of the distribution grid in Section[Tl] The Belief 
Propagation method is detailed in Section [Til] We describe population dynamics algorithms and present SAT-UNSAT transition 
results in IV Control algorithms are explained in Section [V] Section VI summarizes our approach and proposes extensions of 



the study for future explorations relevant to intelligent optimization and control of the power grid. 

II. GRID MODEL 

Consider M sources/generators each connected to D distinct consumers, so that the total number of consumers is MD. 
Greek/Latin indices will be reserved for generators/consumers. For simplicity we assume that each generator has a maximum 
production rate of unity, i.e. y — (0 < y a < 1 jot = 1, • • • ,M), though inhomogeneities in the production can be easily incorpo- 
rated into our approach. The configuration of loads, x = [xi > 0\i — 1, • • ■ ,iV), is drawn from an assumed known distribution 
with support on the interval (0, 1). Our enabling example is the flat-box ensemble with a fraction e of customers dropped off 
from the network, i.e. drawing no electricity at all, 

T( X )=i\[{i-B)p( Xi )+B8( Xi )], p{ X ) = { l/ £ (i) 

The mean consumption then is (1 — e)x, and A is the width of the part of the distribution, correspondent to nonzero demand. We 
require that the consumption of each individual customer is non-negative 1 , i.e. 

2x > A, (2) 

and we assume that under normal conditions, i.e. when demand is not excessive, there is always enough power produced in the 
network to supply electricity to all consumers, i.e. 

M 1 

In our "standard" model of the grid each consumer is assigned to strictly one generator, while each generator feeds exactly D 
consumers. This corresponds to a graph decomposed into M simple trees, where M is the number of generators. To guarantee in 
this standard case that a generator is always SAT, i.e. it is capable of supplying electricity to all of its customers, one requires 

d(x+*)<1. (4) 



1 This assumption is not crucial and can be easily removed, allowing consumers not only to consume but also to generate electricity. The inclusion of this 
capability potentially has a far reaching consequences because it allows distributed generation. We postpone our discussion of this interesting possibility to 
future publications. 
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R = 3 



(a) R = 0, 1,2,3. Graph samples. Ancillary connections to (b) R = 1. Three valid (SAT) configurations (shown in black, 
foreign generators/consumers are shown in color. the rest is in gray) for a sample graph shown in Fig.[TJi. 

FIG. 1 : (Color online) Illustration of the bipartite graph construction of our grid with D = 3. 

This later condition will be challenged below, in the sense that we will show that by adding ancillary lines one does not need to 
impose Q, but instead have a weaker condition, thus extending the domain of the asymptotically failure-free generation. 

To achieve this improvement we consider an interconnected grid built from the standard/separated one by adding ancillary 
lines between consumers and generators. For our quantitative analysis we choose a random graph drawn in the following way. 
First, D consumers are connected to every producer. Second, we choose R consumers from each producer and connect them to 
a second producer in such a way that every producer is connected to D + R consumers. Of these D + R connections D — R go to 
singly connected consumers and 2R to doubly connected consumers. So that a total of NR/D consumers are connected to two 
generators and N(l —R/D) are connected to one generator. See Fig.[T|for illustration. The choice of the two parameter (R,D) 
graph ensemble is used for simplicity. Our analysis method is valid for any random network as long as it is bipartite and locally 
tree-like, meaning that the length of the shortest loop going through a random node is 0(logN). 

We say that given configuration of loads x — [xi\i — 1,* • • ,N) is SAT (satisfiable) if there exists a matching <x — (0;a — 
0, l|{/,a} € Q), where Q is the bi-partite graph accounting for all generators, consumers and lines, such that the following set 
of conditions are simultaneously satisfied 

V/€£: J> te =l, (5) 

Va € Q : £ o ia Xi < 1 , (6) 

ieda 

where di, resp. 8a, stand for all the nodes to which i, resp. a, is connected. If the reverse is true, i.e. there exists no valid <x 
with all Eqs. \5\6\ satisfied, we say that x is UNSAT (unsatisfiable). First, we aim to solve the decision problem: is x is SAT or 
UNSAT? Furthermore, if x is SAT, we would like to find at least one valid solution, <x. Both problems can be stated for a given 
graph or, alternatively, can be considered "in average" for ensemble of graphs. We describe the "average" solution of the former 
problem in Sections IV and the algorithmic solution of the latter in Section [V] Both Sections will be preceded by a preparatory 



discussion of the BP approach in Section[TJI] 

Related models have been studied in the context of resource allocation. One example is the problem of online advertising 
known as AdWords, in particular its uniform off-line version |fl9l l20l |2D . The main difference between our model and the 
model for budget-constrained advertising problem is that condition |6]l is replaced by the maximization of the total revenue 
£ a min[l, Eieto OfocXi] . This relaxation translates into allowing one to shed loads, that is not an option in our setting. Discussion 
of the BP-based approaches to the AdWords problem and their methodology in [22| is very much related to ours. 

Of other models developed in Computer Sciences, the one most similar to ours is the so-called off-line weighted balls-into-bins 
games (sometimes called balanced loading) [23, 24 1. The off-line version of the balls-into-bins games is known to be polynomial 
in the non-weighted case (i.e. for a special discrete choice of the distribution ([T}), as in that case it can be mapped into a max- 
flow (min-cut) problem 11251 . On the other hand the general (weighted) version on a fully connected graph (every customer being 
possibly connected to every generator) is equivalent to the NP-complete bin packing problem (number partitioning in the case 
of two generators). We are not aware of any computational complexity results for the case of bounded number of consumers per 
generator. Therefore, and given that the max-flow mapping of [25] does not generalize to the most general version of Eq. ([T}, 
we conjecture that computationally the general problem described by Eq. ([T]) is NP-hard. 
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III. BELIEF PROPAGATION 

In the asymptotic limit of an infinite system, for which N — > °° while D and R are 0(1), the interconnected grid is locally 
tree-like. Therefore, the Bethe-Peierls (Belief Propagation) BP approach for evaluating the generalization of Q is expected to 
be asymptotically exact (in a sense to be clarified later in this text). This Section describes details of the BP approach 

We introduce the following set of marginal probabilities: 

• v)/™ - *' is the probability that generator a is satisfied given that the edge (z, a), connecting a with his consumer-neighbor on 
Q, is in the active state, i.e. G ia = 1 . 



\|/q^ ( is the probability that generator a is satisfied given that the edge (/,a), where (i,a) £ £71, is inactive, i.e. 0, a = 0. 



• X'{^ a is the probability that i is satisfied (i.e. it is connected to exactly one generator) given that the edge (z, a), where 
(i,a) £ £71, is active, i.e. o, a = 1. 

• Xo^ a i s m e probability that i is satisfied (i.e. it is connected to exactly one generator) given that the edge (/, a), where 
(i,a) £ Q\, is inactive, i.e. o !tt = 0. 

BP relates these marginal probabilities to each other assuming that the relations are graph local, i.e. as if the graph would contain 
no loops. The resulting Belief Propagation (BP) equations are 

xV a = ^ n yf, (7) 

*r = ^ 1 v'r n vr*. (8) 

Pe3;\a ye3r'\a,p 

cr 3c( \ ic( ={0,l}l 3 a\''l /e3a\i ;e3a\i 

yT'' = ^7 E 6(1- E W) IT xC> (10) 

CT 3a\ia={04} |aaV1 jeda\i jeda\i 

where di \ a is a standard notation for the set of generator nodes linked to consumer 2, however excludi ng g enerator a, and 
similarly 3a \ ; stands for the set of consumer nodes linked to generator a excluding consumer i. In Eqs. ( |7|10| , Z'^ a and Z a ^' 



are normalizations ensuring that, %'[^ a + XcP" = 1 an d Vj 1 ^' + Vo^' = 1 ■ The 9( • ) is the step function enforcing the generation 
constraints. It is unity if the argument is positive and zero otherwise. The probability for the link/edge (z,a) to be active, stated 
in terms of the related \|/ and %, is 



p(i,CL) = „_>,• i-y a , . lf a->i v i-»a ' (H) 
Yl Xl +Y0 Xo 

where Y*aedi p(h ot) = 1. Note, for the sake of completeness, that the BP equations can be derived in the spirit of j26l as conditions 
for a fixed-point (minimum) of a functional of marginal probabilities (beliefs), the so-called Bethe free energy functional. 

The Bethe entropy, defined as the logarithm of the number of possible SAT configurations (and also equal to the Bethe free 
energy evaluated at the solution of the previously mentioned BP equations), is 

Sfiethe = Ll0g(Z°) + £log(Z i ) - £ log(Z ! «) , (12) 

a i (j,a) 

z a = 1 e(i-£a ia *,)n&r> (13) 

z' = £ ¥ r ; n t^. (1 4 ) 

aedi |3e3i\a 

Z' a = ^ , zr a +¥o"'Zo" a - (15) 

The entropy Seethe is extensive, 0(N), and self-averaging, i.e. the distribution of Seethe is concentrated around its mean value 
with dispersion being o(N) at N — > 00. 
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IV. AVERAGE BETHE ENTROPY VIA POPULATION DYNAMICS 



Fixed point of the BP Eqs. (TJlOl and the corresponding entropy ( 15 1 can be obtained by solving Eqs. (7|[T0i iteratively for a 



given instance of the problem. Repeating these simulations many times for different instances of the graph and load ensembles 
one can also calculate the average Bethe entropy. However, the average behavior, corresponding to the limits of the infinite 
graph, can be obtained more efficiently via the population dynamics technique. This technique offers a computationally efficient 
sampling from the distribution of the marginal probabilities % a on infinite random graphs and subsequent evaluation of the 
average Bethe entropy. Below we give a very brief exposition of the population dynamics approach applied to our model. The 
interested reader is referred to [ 1 1 , 27 1 for further details. 

In population dynamics, we create a pool of N p components, each characterized by the (5Co,5Ci,*) vector with x drawn from 
the original distribution of demands ([T]) and the initial % selected arbitrarily. Messages leaving the generator and entering the 
consumer connected to a single generator will always be kept fixed to Xi = 1, Xo = 0, while other messages will be updated 
iteratively, so that at any new sweep a new pool is derived from the old one. A sweep consists of the following step repeated 
0(N P ) times (each time one of the N p components is updated). A step consists of choosing a random number x,, D — R random 
numbers xj, an d 2R random elements from the pool representing the incoming messages and computing the message \|/ based 
on Eqs. ( 9|10 l. Then we replace a random element of the pool by the vector ,\|/o,x,), in accordance with Eqs. (7J8 I, thus 



guaranteeing that %\ = \\Iq and %o = We repeat this sweep procedure many times until convergence is achieved. To compute 
the average Bethe entropy we apply similar procedure of sampling from the resulting pool and thus averaging all the terms in 
Kq. {15]). 

Implementing the population dynamics method we observed three possible outcomes 

(a) SAT phase: The Bethe entropy is positive suggesting that the number of SAT configuration (valid redistribution of the 
demand over the generators) is exponential in the system size 2 . 

(b) UNS AT phase, type 1 : The Bethe entropy is negative, suggesting that there is almost surely no valid redistribution of 
demand over the generators. 

(c) UNSAT phase, type 2: A contradiction is encountered in the BP equations, formally correspondent to zero values for 
the normalizations in (TJlOl. We conclude that the demand is incompatible with the graph and respective generator 



assignment. 

Fig. [2] shows results of the population dynamics simulations, with lines connecting the marks indicating the boundaries of the 
respective SAT and UNSAT phases. These boundaries are established by fixing A and traversing different values of x starting 
from SAT phase moving towards UNSAT phase and catching the value where the UNSAT conditions are first time observed. We 
observe a significant improvement of the SAT/UNSAT threshold from what Eq. Q suggests. 

Recall that in the R = case the SAT/UNSAT threshold is described by fi). Similar upper bound can be derived for cases 
R > 1. Then for any given R, one can find (with high probability) a place on a very large network where all consumers connected 
to two generators (which are neighbors trough one of the consumers) have demands near to the maximum value x + A/2. Then 
at least D — R+l of these loads have to be connected to one generator and hence 

(D-R+l){x+~)<l. (16) 

Note that ([16]) is identical to Q for R = 1 . 

Moreover, to observe an improvement in the condition Q at e = one needs to require existence of such configuration of 
loads that a set of (D+ 1) consumers connected to one producer and drawing the minimal amount of electricity does not overload 
the generator, i.e. 

e = 0: (D+l)(x-~)< 1. (17) 

To conclude, this study of the average Bethe entropy shows that the network with added ancillary lines is able to withstand 
larger fluctuations in the demand, A, than the naive network, in which every consumer has a pre-designed provider independently 
of the current demand. This effect is amplified with increasing R. 

At this point it is also appropriate to recall that the SAT-UNSAT transition is actually an abrupt transition only in the sense of 
the asymptotic N — > °° limit. Thus for large finite N the generator failure probability is small but finite at any point of the SAT 
domain. 



2 This is the case unless the so called replica symmetry breaking takes place in our model. We have done a local stability check of the BP solution and have not 
seen any indication for a break down of the replica symmetry. Note also that this stability is also equivalent to convergence of the underlying BP algorithm 
for an individual realization of the graph and the loads. 
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FIG. 2: (Color online) Results of the population dynamics shown in the (x,A) plane for D = 3,4, where (1 — e)x is the distribution mean 
and A is the distribution widths. The population dynamics was done with Np = 1000 for D — 4 and N p = 10000 for D = 3. The triangular 
region (black full lines) correspond to conditions {2|, {3} and (T7}, and encloses the interesting region of parameters. The colored lines with 
data points separate SAT and UNSAT domains lying on the bottom-left and top-right off the lines respectively. Different colors/lines/markers 
correspond to different values of R = 0, • • • , D. Note that the curves for R = and R = 1 (dashed without markers) are identical, due to condition 
\16) . We observe that the performance improves with increasing R>2 (the number of ancillary lines) and the grid tolerates larger values of A 
(fluctuations in the demand). The abrupt change in slope of the curve corresponding to R = 2 (full red line with x marks) in the right part of the 
figure is due to condition l |16| >. Note that if condition |2j is removed (see e.g. discussion in the footnote preceding the equation) our description 
still remains valid, thus resulting in the colored lines with markers extending smoothly beyond their crossings with the tilted black curve. 



V. CONTROL ALGORITHMS 



In this Section we discuss the problem of calculating a valid load-to-generators assignment for a given graph and given 
configuration of loads. We designed two heuristic methods to identify SAT configuration of switches. Our first algorithm, coined 
WalkGrid, is an adaptation of the WalkSAT l28l|29l , which is a stochastic local search heuristic solver for the K-satisfiability 
problem. The WalkGrid algorithm is very fast and shows flawless performance in discovering a valid configuration almost 
anywhere in the SAT region. 

Our second algorithm corresponds to solving BP equations. The BP scheme, designed in the spirit of fl2l . is used to find 
most biased/stressed link and then proceeds with decimation towards a valid configuration [12|. Our implementation of the 
BP-decimation is so far slower and a bit less efficient than performance shown by the WalkGrid algorithm. To this point, let us 
note that there exists a more efficient way of using BP to find valid configurations - the reinforcement strategy [30 , 31 1, which 
is fully distributed, linear in the number of consumers and typically outperforms decimation. The reinforcement strategy has 
been implemented in ll22l for a related online advertising problem. However, we had a difficulty to find an implementation of 
the reinforcement which would work efficiently in our problem. 



A. WalkGrid 

WalkGrid is a stochastic local search algorithm inspired and closely related to its K-SAT ancestor called WalkSAT Il28ll29l . 
It can also be viewed as a Monte-Carlo-like algorithm which in order to gain speed violates the detailed balance condition (not 
needed here as we are not interested in sampling, but are rather focused on a local search). Our implementation works as follows 
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WalkGrid 

1 Assign each value of o or 1 randomly (but such that V/ 6 G : Y.aedi'^ia = 1); 

2 repeat Pick a random power generator a which shows an overload, and denote the value of the overload, 8; 

3 Choose a random consumer i connected to the generator a, i.e. G, a = 1 ; 

4 Pick an arbitrary other generator which is not overloaded and consider switching connection from (/a) to 

5 if (in the result of this switch a is relieved from being overloaded 

6 and p either remains under the allowed load or it is overloaded but by the amount less than 8) 

7 Accept the move, i.e. disconnect i from a and connect it to p thus setting o,p = 1, G ltx = 0. 

8 else With probability p connect consumer i to p instead of a; 

9 until Solution found or number of iterations exceeds MT max . 

The WalkGrid algorithm depends on two parameters: the maximum number of iterations T max , and the temperature-like 
(greediness) parameter p. The parameter p needs to be optimized, just as in the original WalkSAT solver. 
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FIG. 3: (Color online) Performance of the WalkGrid algorithm. The data are for networks with M producers, and R = 2, D = 3, A = 0.2, e = 0. 
The temperature-like parameter and the maximum number of iterations are set to p = 0. 18 and T max = 2000 respectively. The average (median) 
is over 100 instances for M up to 20k, and 20 instances for M = lOOfc. Left: The median running time plotted against the average consumption 
x. Right: The percentage of cases where solution was found in less than T m . dx = 2000 iteration. The Bethe entropy based (asymptotic) curve 
is drawn dashed for comparison (the actual value for the curve is not related to the success rate), suggesting that in the limit of N — > °° valid 
configurations exist up to x ~ 0.30. Note that in the separated case (of R = 0), valid configurations exist only atx< l/D — A/2 ~ 0.23. 

Fig.[3]shows performance of the WalkGrid algorithm. The running time of this algorithm scales close to linear with the system 
size, it is thus relatively easy to resolve fast network with many thousands of nodes. Whereas for R = 2, D = 3, A = 0.2 the 
separated architecture requires x < 0.233, the WalkGrid algorithm is able to find valid configurations up to x ~ 0.296, while our 
theoretical analysis suggests that valid configurations should exist up to x « 0.301. 



B. Belief-Propagation decimation 

In the BP-based decimation algorithm one updates Eqs. (|7]|T0]> iteratively, thus passing messages from generators to consumers 
and back. After fixed number of steps the most biased consumer is chosen and the moreprobable value for its consumption is 
assigned, the graph is reduced and the procedure is repeated. Note that updating Eqs. (|V]|8j» takes 2 2R steps per message thus 
making the algorithm exponential in R. However, building new connections is expensive and one should realistically assume 
that actual R is any case not very large. As far as the scaling in goes, the algorithm is quadratic in the number of consumers. 

Decimation 



1 repeat Update BP messages on every edge according to (7fll0i n times 



2 Compute the marginals (Hi 

3 Choose the most biased edge ia and assign is the more probable values; 

4 Simplify the formula, cutting off the assigned edge from the graph; 

5 until Solution or contradiction is found; 
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FIG. 4: (Color online) Performance of the WalkGrid algorithm for D = 3, R = 2, e = and different values of the distribution width A. Bellow 
the blue (dashed) line the WalkGrid algorithm is able to find solution on more than 50% (average over 20 trials) of graphs with 10 5 generators. 
The red (full) line corresponds to the theoretical boundary separating SAT and UNSAT domains respectively. 



The algorithm performance is illustrated in Fig.[5]where the percentage of success in the BP-based decimation is shown. These 
data average over 50 random instances from the R = 2, D = 3 ensemble with variance in consumption A = 0.2 and varying mean 
consumption. For such a set of parameters the separated network would be able to achieve x < 1/3—0.1, whereas with two 
consumers per producer connected to two producers the mean consumption increases to x ss 0.30. 




0.32 
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FIG. 5: (Color online) Performance of the BP-based decimation algorithm. The data are for networks with M producers, n = 5 (five iteration 
per a cycle of the decimation procedure) and R = 2, D — 3, A = 0.2, e = 0. Average over 50 random instances is taken and the fraction of 
successful runs is plotted against the mean consumption x. The Bethe entropy based (asymptotic) curve is drawn dashed for comparison (the 
actual value for the curve is not related to the success rate), suggesting that in the limit of N — > 00 valid configurations exist up to x ~ 0.30. 
Note that in the separated case (of R = 0), valid configurations exist only at x < l/D — A/2 ~ 0.23. 



VI. SUMMARY AND PATH FORWARD 



This manuscript reports a first study of the power distribution networks with transition from the SAT regime, in which shedding 
of loads is avoidable, to the UNSAT regime, where shedding is the only available option for balancing the demand. We have 
shown that a significant enlargement of the SAT-domain is possible employing ancillary connections between power consumers 
and generators. Even though our model represents a gross oversimplification over the actual power grid, it offers a significant 
step forward in providing a framework and guidance for analysis of more involved and realistic problems. The general approach 
we pursued in this study is based on recent developments in the field of graphical models that merges statistical physics, computer 
science, optimization theory and information theory lfl3llT41 . The BP approach is asymptotically exact on infinite sparse graphs 
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and as such is useful for the asymptotic (capacity/phase-transition style) analysis. The BP scheme also provides heuristic tools 
for graphical models on finite sparse graphs that can be used for algorithmic optimization and control of the power grid. We 
tested this BP approach, and also developed in parallel another an apparently more efficient alternative to BP called WalkGrid. 
This algorithm finds valid configuration of switching practically anywhere inside the SAT phase. Note, that the two algorithms, 
BP and WalkGrid, are truly complementary and there utility for practical problems in power networks are yet to be explored. 

It is important to emphasize that many generalizations of our model are very straightforward and can be used directly within 
the framework presented here. This includes implementing different probability distributions for demands (as long as the dis- 
tribution support is bounded), and non-uniformity for both generator and consumer levels, i.e. varying production caps for 
generators and introducing distinct distributions of demands for different consumers. Also the network itself can be easily 
extended in many ways from the simplified case, parameterized by R and D, that we discussed in the manuscript. Our equa- 
tions are straightforwardly valid for every bipartite locally tree-like random network. Another case, which allows very natural 
and straightforward generalization for all the statements made in this manuscript, corresponds to breaking the equivalence be- 
tween different edges in the graph and thus assigning nonuniform weights to them. These weights may e.g. represent cost of 
construction, geographical length, proxy for losses, cost of exploitation, etc. 

There are many other more realistic extensions of our model associated with the description, optimization and control of 
power grids which can benefit from utilizing a graphical model approach of the kind discussed in this paper, even though actual 
implementation may prove to be more involved. We conclude listing some of these more interesting but difficult problems that 
we plan to address in the future, based on the general method sketched in this manuscript: 

A. Most important generalization of our approach would be to account for losses, impedances and the reactive character of 
AC electrical systems. Obviously this will require incorporating in our statistical SAT-UNSAT framework Kirchhoff's 
circuit laws Q [32] . 

B. Generators in a real grid are interconnected on higher (still power distribution, but also power generation) levels. In com- 
bination with item [A], this represents a major challenge for extending our approach. However, we still believe that posing 
the joint optimization and control problem in terms of a complex graphical model and then addressing capacity /transition 
as well as algorithmic issues with the host of BP-related techniques is a feasible and exciting path forward. 

C. The joint optimization setting, mentioned in item [B], may also include various additional factors associated with economic 
policies (e.g. prices, incentives, etc) Il9l l33l . government regulations and load/generation forecasting. These complications 
can be accounted for in the form of extra soft or hard constraints. This united framework should also take advantage of 
the progress made in developing optimal power flow solutions |j34ll35l . 

D. Our consideration in this paper was purely static, thus ignoring important transients. Constructing dynamical (discrete 
and continuous time) models that also accounts for all the aforementioned problems, is an important future task. It is also 
important to notice that the graphical model approach can in fact be extended to the dynamic framework, see e.g. [36|, 
thus suggesting yet another intriguing future opportunity. 

E. An essential part of statistical studies of the power grid focuses on estimating probabilities and mitigating very costly and 
dangerous large scale outages ll37l [38l . The graphical model framework this manuscript describes allows an extension 
which can analyze rare events of special interest, such as dangerous but rare configurations of the channel noise and 
degrading performance of Low -Density-Parity-Check codes [39 40]. We plan to extend this approach to the analysis of 
outages and their cascading through the grid. 
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